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ABSTRACT. Newcastle disease is common viral poultry disease which leads to a massive 
killing of chicken if preventive measures are not well taken. In this paper, we develop and 
analyze a deterministic model to investigate the role of wild birds and environment on 
the transmission dynamics of Newcastle disease in village chicken. We compute the basic 
reproduction number (Ro), a threshold that tells the presence of the disease in a population. 
Finally, we performed the sensitivity analysis of parameters to see their relationship with the 
basic reproduction number (Ro). The numerical results show that the basic reproduction 
number (Ro) is more sensitive to the contact rate between the susceptible village chicken, 
wild birds and contaminated environment. This implies that, more contamination of the 
Newcastle virus into the environment increase the chance for the repeatedly occurrence of 
the disease. The results also shows that increasing the clearance rate of Newcastle disease 
virus in the environment reduces the rate of spread of the disease in chicken population. 
Therefore, contaminated environment plays a crucial role in the transmission of Newcastle 


diseases in the village chicken population. 


1. INTRODUCTION 


Newcastle disease (ND) is a viral disease that affects poultry and other wild birds [3, 10]. 
This disease is a major constraint to the development of village chicken industry particularly 
in Africa and Asia since it causes a high morbidity and mortality rate of chicken [3, 13, 24, 28]. 
The Newcastle disease is caused by Avian paramyxomovirus-1 (APMV serotype-1) virus in 
paramoxyviridae family |13, 22, 27]. Depending on the severity of the disease during the 
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outbreak, NDV is characterized as lentogenic, mesogenic or velogenic strains [4, 27]. Velo- 
genic strains cause acute disease with high mortality of up to 10096 for un-vaccinated village 
chicken compared to other pathotypes [27]. 

The symptoms of ND includes paralyzed legs and wing of the affected avian, coughing, head 
twisting, watery and greenish diarrhea and other nervous symptoms that follow in one or 
two weeks since infections [2, 25]. T'he seriousness of the disease depends on the virulence of 
the infecting viruses, susceptibility of the host, other diseases in the flock, environmental in- 
fluences, and the vaccination history of the chicken [3]. Young birds in a flock are extremely 
susceptible to disease where death rate is about 100% [16, 24]. 

Although the Newcastle disease is endemic and sometime epizootic in village poultry, many 
aspects of its epidemiology are not clearly understood[13]. According to [19], other domes- 
tic birds, carrier chicken, wild birds and the physical environment are thought to be the 
contributing factors in the cycling occurrence of Newcastle disease. Initially, the disease is 
transmitted when an infected chicken is introduced in the flock. Further contacts of the 
infected chicken with other un-affected chicken leads to the spread of the virus in the entire 
flock [2, 27, 10]. Inhalation and ingestion of the virus through water and feed are the primary 
means of the spread of the virus [24, 27]. The infected and carrier birds may also shade NDV 
in their discharges and contaminate the environment. Depending on the season, the virus 
can survive for days in the environment, forage, water and in the birds feathers [4]. The 
contact of the susceptible birds with the contaminated water, food, droppings or discharges 
of the infected birds and other farm utensils can also lead to the spread of the virus [8]. 
Over the years, different scholars have formulated and developed mathematical models to 
gain insight on the transmission dynamics of Newcastle disease in the village chicken|20, 14, 
6]. 

[20] extended the model for foot and mouth disease in Thailand [29] to model the dynamics 
of ND among the village chicken. The model assumed that village poultry are the reservoir 
of Newcastle disease virus for other poultry sector. The model revealed that, in order for the 
vaccination of ND to be effective it must be conducted frequently in a large population of 
the village chicken. [14] developed a deterministic compartmental eco-epidemiological model 
with optimal control of Newcastle disease (ND) in Tanzania. The study revealed that ND 
can be controlled if farmers vaccinate chicken properly in time. [6] developed two mathe- 
matical models to study the influence of illegal harvest and effects of age structure of the of 
wild wingled-Parakeets on the dynamics of Newcastle Disease. The relationship between of 
Newcastle Disease transmission and harvest were evaluated through the basic reproductive 
number. The findings showed that much harvest reduce the spread of the ND in the Para- 
keets population. Using their second model they concluded that age of the Parakeets can 
influence the spread of ND though not to a great extent. 

To the best of our knowledge no any mathematical model which tried to address the role 
of wild birds and environment in the transmission of Newcastle disease. Therefore, based 
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on different clinical and epidemiological studies conducted on the transmission of Newcastle 
disease in village chicken, we develop a mathematical model to investigate the role of the 
environment and wild birds in the transmission dynamics of Newcastle disease in village 
chicken. 


2. MATERIAL AND METHODS 


2.1. Formulation of the Model. We develop a compartmental model where chicken and 
wild birds are categorized in different compartments depending on the status of their infec- 
tion. The Chicken population are modeled as S,E.I. (Susceptible-Exposed-Infected) while 
wild bird population are modeled as Sy EjI,I, (Susceptible-Exposed-Infected-Carrier). The 
environment has only one class of Newcastle disease virus (NDV). The model is formulat- 
ed based on the idea that wild birds and the environment are the primary reservoirs of 
Newcastle disease virus [2, 10, 24, 19, 17]. The village chicken population N,(t) is divided 
into three subpopulations namely: S.(t) denoting the number of susceptible village chicken, 
E(t) denoting the number of the exposed village chicken, and /.(t) denoting the number of 
infected village chicken. The Wild birds population Nj; (t) is divided into four compartments 
namely; S (t) denoting the number of the susceptible wild birds, Ej(t) denoting number of 
wild birds exposed to NDV, I(t) denoting the number of chronically infected wild birds and 
[,(t) denoting the number of reservoir wild birds. The environment contaminated with NDV 
is denoted by H(t). In the model, we assume a density dependent recruitment rate uN, for 
the village chicken population through birth. Initially, village chicken acquires Newcastle 
disease virus when a diseased village chicken is introduced in a flock and come into contac- 
t with other un-affected chicken. Village chicken can also acquire Newcastle disease virus 
when exposed to unhygienic environment as well as when interacting with other carrier wild 
birds [17, 6]. Thus chicken acquire infections at a rate qi(t, Ie, L,, H) that depends on per 
capita contact rate of the infected village chicken with its susceptible population (vI.(t)), 
per capita contact rate between the carrier wild birds and the susceptible village chicken 
eu). The infect- 
ed chicken becomes infectious after surviving within the incubation period (1/7) of two to 





(bI,.(t)) and the density of Newcastle disease virus in the environment ( 


six days since infections |2]. After surviving the incubation period, chicken starts to show 
aerosol signs and progress to chronic stage of the disease at a rate yE.(t). We assume that 
village chicken cannot recover from the infections but they die naturally at a rate u and by 
the disease induced death rate ôe. 

Wild birds are assumed to be recruited at a density dependent rate uM, through birth. Like 
the village chicken, the susceptible wild birds gets NDV from the contaminated environment 
as well as when they interact with the clinically infected and reservoir wild bird population 
at the rate y(t, Ia, L., H) which is a function of the interaction between the susceptible wild 
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birds and the clinically infected wild birds (¢/,(t)), the reservoir population of wild bird- 
s (al,(t)) and (ae), where the parameter d is the contact rate of the hosts with the 
environment. The latently infected population of wild birds is considered as a reservoir of 
the virus[4]. This makes the infected population of wild birds to be divided into two sub- 
populations; chronically infected, [,(t) and carrier, [,(t),at the proportions of p and (1 — p) 
respectively. It is assumed that wild birds cannot recover from the disease once infected but 
they are reduced by the natural death at the rate u and by the disease induced death at the 
rate 0y.The Carrier wild bird population does have disease induced mortality, its assumed 
that they only die naturally. 

Newcastle disease virus is introduced in the environment through shedding by the clinically 
infected village chicken, the clinically infected wild birds and the carrier wild birds at the 
rates a, and ay» respectively[4, 24] . The Newcastle disease virus can survival for some months 
at a temperature between 20 — 30°C and much longer at cooler physical environments [19]. 
Parameters and state variables used in the formulation of the Newcastle disease model are 
summarized in the tables below; 


TABLE 1. Descriptions of parameters used in the formulation of the model 
system 2.1 





Parameter Description 


Selt) Susceptible village chicken population 
E(t) Exposed village chicken population 

TX) Infected village chicken population 

Salt) Susceptible wild bird population 

E(t) Exposed wild bird population 

I(t Chronically Infected wild birds population 


Carrier wild bird population 
t) Total population of village chicken 
t) Total population of wild birds 
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TABLE 2. Descriptions of parameters used in the formulation of the model 
system 2.1 





Parameter Description 





a Transmission rate between carrier wild birds and the susceptible 
population of wild birds 

b Transmission rate between carrier wild birds and the susceptible village 
chicken 

Ot Shading rate of NDV in the environment by chronically infected 


and the carrier wild birds 
Qe Shading rate of NDV in the environment by the chronically infected 
village chicken 


T Recruitment rate and Natural mortality death of the host populations 

[o Transmission rate between chronically infected wild birds and the 
susceptible wild birds 

v Transmission rate between chronically infected village chicken and the 
susceptible chicken 

K Half saturation constant of Newcastle disease virus in the environment 

d The contact rate between the susceptible population of village chicken 
and wild birds with contaminated environment 

p Proportion of the exposed wild birds which become chronically infected 
with NDV 

hia Clearance rate of the NDV from the environment 

Ob Disease induced death rate in wild birds 


ĝe Disease induced death rate in the village chicken 
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2.2. Model Flow diagram. Based on the transmission dynamics of the Newcastle disease, 
model assumptions, definition of variables and parameters respectively, the dynamics of the 
Newcastle disease is summarized in the flow diagram as follows: 





FIGURE 1. Compartment model diagram for the transmission dynamics of 
Newcastle disease in village chicken population. The sold lines shows the con- 
stant transmission from one compartment to another, the dotted lines shows 
the normal interactions between different compartments and a dash dot lines 
represents the shedding of NDV onto the environment 
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2.3. Model Equations. The dynamics of NDV is described by the following system of 


nonlinear differential equations; 


























Ll. = uNe — p1 (t, Ie, L,, H)S.(t) — uS. (t), 
= = eit; Ie, Ir, H)Se(t) — (y + p) Eelt), 
ze = yE,(t) — (ôe + n). (t), 
lO = uN, — polt, Dp, Ip, H)S,(t) — S(t), 
(2.1) po zs pae d RE RSS CO Sein Bey. 
UU = py Ey(t) — (0s + w)Lo(t), 
d = (1 —p)yE,(t) — uL(t), 
UD = palt, Io I, Tr) — to H(t) 
where; 
eit Ie Ie, H) = (52 yt aoe z i) 
óLy(t) + al,(t) H(t) 


pat, Ij Ip, H) = ( 








No 


p3(t, d l, 1,) E a l.(t) + Ob (L(t) F [,(t)) 


With initial conditions; 


K+ A(t 


jJ 


S.(0) > 0, E.(0) > 0, 1.(0) > 0, S,(0) > 0, E(0) > 0, I(0) > 0, L,(0) > 0, H(0) > 0. 


2.4. Basic properties of the model. 


2.4.1. Invariant region of the solution. The Newcastle disease model system (2.1) has three 


subpopulations where all parameters and variables are positive Vt > 0. 


Lemma 2.1. Given the model system (2.1) in 














RÈ with the initial conditions S.(0) > 0, 


E.(0) > 0,1,.(0) > 0, S,(0) > 0, £,(0) > 0, (0) > 0, (0) > 0, H(0) > 0, its solution 
enters the invariant region 


Q = (S«(t), E(t), I.(t) € 


























R3 Selt), Ext), I(t), I(t) € 


R4; H(t) € 














uj) 


Proof; To establish the feasible region of the Newcastle disease model solution, we apply the 














box invariant method as used by [1, 15]. For our dynamical system g = g (X,t), X € R”, we 
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assume the continuity and the Lipschitz properties of its solution. The model system (2.1) 
is reduced to the form 


dX 
2.2 — = XG 
(2.2) <= Q(a)X + 
where X = (Se, Ee, Io, S Eo, Io, do H)” and a column vector G = (Na, 0,0, Ny, 0,0,0,0)” > 
0. G is a Lipschitz continuous[1, 21] and 


(2.3) Q(z) = 0 Q(z) 0 
0 0 Q3(z) 














is a Metzler matrix for all X € RẸ with sub-matrix Qi (x) , Q» (x) and Qs (x) from the 
village chicken, wild birds and environment respectively. We define the sub matrices from 
the system (2.3) as follows; 


— (gilt, Lo Ir, H) + u) 0 s 
(2.4) Qi (x) = yi(t, Ie I, H) "ite 0 
— (volt, I5, I., H) +u) 0 0 X 
os q»(t, I», fp, H) ~(y +1) 0 ki 
(2.5) Q2 (x) i " is A F ; 
0 (1 a p)y 0 —u 
(2.6) Qux) - (0 0 as 0 0 ap a n.) 


By combining the matrices in equation (2.4) , (2.5) and (2.6) we get the matrix Q (x) which 
is a Metzler matrix for all X € R$ defined as; 














-A 0 x 0 0 0 0 0 
A, —(w- u) 0 0 0 0 0 0 
0 y (b+) 0 0 0 0 0 
b 
ena@=|, 4 ES a 4E M 
a —(¥+H) 0 N; 0 
0 0 0 0 py —(óp du) 0 0 
0 0 0 0 G-p 0 —u 0 
0 0 a 0 0 a Q  —pwK 
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where 


gilt, IL, I, H) Fu 
Oil Ie, Ip, H) 
v»(t, D, I, H) + H 
Ag = alt, Ip, Ip, H) 


Ay 
Ay 
As 


A reduced Metzler matrix Q (x) in (2.7) has all negative values along its principle diagonal 
and the rest non-negative values in its off diagonal. Hence proves that all variables enters 
and remains in the invariant region Q. This shows that the Newcastle disease model system 
(2.1) is epidemiologically meaningful and well posed in 2. 


2.4.2. Positivity of the solution. 


Theorem 2.1. Let the initial value of variables of the model in the equation (2.1) be S.(0) > 
0, E.(0) > 0, I,(0) > 0, S,(0) > 0, E;(0) > 0, 7,(0) > 0, 1,(0) > 0 and h(0) > 0 

then the solution set 

{(S-(t) > 0, Eelt) > 0, I(t) > 0, S (t) > 0, E (t) > 0, I, (t) > 0, I(t) > 0, H(t) > 0) € RE is 
positive for all time t. 























Proof 
Lets consider the first equation of the model (2.1) 
dS. (t 
(2.8) D ) = uN. — gilt, Ie, Ip, H)S«(t) — Selt) 
dS.(t) 
29) AD) > — (galt LL, H) + 1) Slt) 
t dS. (t) t 
(2.10) EE (¢1(t, Le, Ip, H) + u) dt 
(2.11) St) 2: S(0)e (foe eterna 


Thus as t — oo then it follows that Se(t) > S,(0)e-^t-JoGn(11-.H > 0, From the second 
equation of the model (2.1) we have 








(2.12) wat = qt, Ie, Ip, H)Se(t) — (y+ u)t 
(2.13) DE) > —(y* p)E«(t) 
(2.14) dES( 5s _ Gy + u)dt 





E(t) 
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Integrating both sides of equation(2.14) w.r.t we have 


(2.15) f n 2 = (y + u)dt 





and finally we get 
(2.16) E(t) > E.(0)e-6*?*t 


As t — oo, Eelt) > E.(0)e-^*?* > 0 , we have E,(t) > 0 
Also from the third equation of the model (2.1) we have 








(2.17) Ml) — YEE) — (6. + HL) 
(2.18) i -u > - fice mat 


which gives I,(t) > I.(0)e~ +")! > 0. Following the same procedures it follows that; S(t) > 
0, E(t) > 0, L(t) > 0, I(t) > 0 and H(t) > 0 which proves that all state variables are 
positive Vf. 


2.5. Existence of the disease free Equilibrium point. 


2.5.1. Disease Free Equilibrium point ( DFEP). 'The model has a disease free equilibrium 
point (DFEP) which is obtained when all infectious states of the model are set to zero. 
Setting /,(t) = I(t) = (t) = H(t) = 0 in the model system (2.1) we then have the disease 
free equilibrium point D? = (N.,0,0), D? = (N,,0,0,0] and D$, = 0 for village chicken,wild 
birds and the concentration of NDV in the the environment respectively. Generally, the 














disease free equilibrium point of a model is given by D? = { N., 0,0, N,,0,0,0,0]. 


2.6. The basic Reproductive number, Ro. The basic reproduction number Ro, is de- 
fined as the average number of secondary cases caused by one infectious individual intro- 
duced in a population that consisting of entirely susceptibles|9, 11, 23]. This number tells 
and quantifies the ability of an infectious disease to invade a purely susceptible population 
and persist|9, 11]. The Epidemic persists when Rp > 1 and dies out when the Ro < 1 
[7, 12, 31, 32]. We compute the Ro by the next generation method as proposed by [31]. By 
[31], we define our system for infectious compartments as 

dX; 
dt 


From equation (2.19), X; defines a set of infected classes, F; defines the rate of new infections 


(2.19) 





=Fi-V; 


in compartment i and the total transfer rate is defined by V; = V; —V;", where V; is the rate 
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of transfer of individuals out of compartment i and V* is the rate of transfer of individuals 
into compartment i through interactions. Then it follows that; 











Ic i Ir | H 
Ue (ve FON 4 dz) Se —(y + i) E. 
de 0 ye E (ôo t He 
“a _| (Se +a) |o —(y+ Wer 
t 0 py Es + (ôb + u) 
di 0 (1 — p)y Ey, 
Hu 0 Qele + Ay (1, + L) — uH 


The corresponding Jacobian matrices of F and V are the matrices of the derivatives of F; 
and V; w.r.t Felt), Ie(t), Ealt), I(t), I(t) and H(t) at the disease free point, D? which are 
given by 


- EOD (VD?) ; 
F= ( aX, ) and V = ( DX, respectively. 








Then we get the following matrices 


0 $ 0 0 bye d^ 
0000 0 0 
(2.20) REX E - 
0000 0 0 
0000 0 0 
0000 0 0 
y+u 0 0 0 0 0 
—y aby 0 0 0 0 
(2.31) Ve 0 0 (y t p) 0 0 0 
0 0 —py (+p) 0 0 
0 0 == pyr 0 u 0 
0 =e 0 —Op —Ay fy 
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The inverse of matrix V in equation (2.21) become 














(yu) 0 0 0 0 
est 
GG) (e) 0 ? : 
0 0 (yu)! 0 0 
(2.22). eom » 4 

0 0 Gun) (t 0 

_ (+p) -1 
0 0 ortu) ? x 
"y Oc Qe _ a(=7 50-7 H+7 p Ôb) a on 
(+u) (de+H) Hv (Oc+H) Ho (+u) Oo+H) HE Hv (do +H) Ho nm 


We then compute the next generation matrix (FV 1) which gives 

















dN a bN, | dNeay, dN 
RS T EOE Now | "um K ey 
0 0 0 0 0 0 
9 , dN a^ rdi ue dN 
(2.23) FV! = Ge Ee dotu | KOED 7 mm Kd 
0 0 0 0 0 0 
0 0 0 0 0 0 
0 0 0 0 0 0 
where 
= vy : AN oy a, 
(y +u) (Setu) r (Y+ u) (6s + B) us 
EN ME dN Qe 
EDEN | K (e + H) by 
po Nell +p)y  dNow(—yÓ, 1H + 70 56) 
No (y+ #) B K (Y + u) (ôe +B) b p 
dN yy Ae 





(V+) (ôe + p) be 


u AN yo. 
K (de + H) My 


So Oo a] C © 
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_ $ py | a(-1-p)w  dNyw(—-yó -YH - vp) 
(+u) (+u) | Qq-unu K (y+ u) (55 + p) H Hov 








Basing on the matrix (2.23), the basic reproductive number Ro is the spectral radius of the 
next generation matrix (p (FV-!)) . This gives the basic reproduction number as 




















(2.24) Ro =p(FV") = ; (œ FW)4 V(R- Wy AUT) 
(2.25) 
= py | a(-1+p)7  dNso(—yôe— yp +Y p) 
Ro =1/2( (wc) -p (Vu) k (y+ u) (o + B) H tw ) 
S py _a(=1+p)7 , dNy(7yà — yn yp) V. | 
ey (n (v-u)(--u)  (w—mnu n (Y+ u) (ô+ pup. ) dt 
for 





E 4dNyy a, (E a sro] 
K (y + u) (de + p) Mo Matau K (y + u) (o + u) H ps 

From the equation (2.25) the basic reproduction number 2 is influenced by parameters 

respectively refers 


from all subpopulations of the model. In this equation, Gu and 


135 13s 
the probabilities that village chicken and wild birds survive throughout eds whole incubation 
period in the presence of Newcastle disease and Tag is the probability of village chicken to 
acquire NDV when come into irae with the icons village chicken throughout its entire 
infectious period. The rate Gn ^ is the probability of wild birds to acquire NDV when come 
into contact with chronically affected wild bird through its entire infectious period. However 
em indicates the probabilities of village chicken and wild birds to acquire the Newcastle 
disease virus from the environment during the outbreak of the disease. 


2.7. Sensitivity analysis of Ro. The sensitivity analysis of y helps to understand the 
behavior of parameters on the model output as well as their influence in the spread of the 
disease in the population [26]. In order to get the sensitivity value of a parameter in the 
model system (2.1) , a function Ro should be differentiable at a parameter q where q being 
any parameter in Ro. Therefore, in this section we perform the sensitivity analysis of Ro 
by using the normalized forward sensitivity analysis index as applied in [5]. A normalized 
forward index of a variable q is defined by 

ORo q 

I Ro 


Parameter values to be used are derived from the field study, the related literatures and 


Ro _ 
ae = 


others are assumed. The population of village chicken in Bagamoyo and Kibaha districts 
of a Coastal region (Tanzania) is estimated to be four hundred thirty one thousand seven 
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hundred and fifty (431, 750) with 90% death cases on absence of any control during the out- 
break of Newcastle disease. This gives an estimate of the Newcastle disease induced death 
ôe = 0.01984day !. Wild birds are assumed to be five hundred thousand (500,000) and its 
disease induced death is assumed to be 6, = 0.0025day~!. We estimate an average lifespan 
of village chicken to be two to ten years which yields a mortality rate u = 1.37 x 10 ?day^! 
[18]. The parameter & gives the density of Newcastle disease virus in the environment that 
yields 50% chance of infection which is assumed to be 10000 virus per cubic meter of the 
soil. Other parameter values and the sensitivity analysis values are summarized in Table 3 


and table 4 respectively. 
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TABLE 3. Parameter values of the model 2.1 








Parameter Parameter value Source 
a 0.01day^! Assumed 
b 0.21day ! [2, 8] 
Me 0.01667virus-!chiken-!day ^! Assumed 
ab 0.02virus lchiken-!day ! Assumed 
V 0.0083day ^! Assumed 
u 1.37 — 5.48 x 10-*day-! [18] 

D 0.02day~! Assumed 
p 0.9 Assumed 
d 0.001day^! Assumed 
y 0.067day~* [29, 30] 
Op 0.025day~* [6] 

à. 0.01989day-' [14] 
Lo 0.00219day Assumed 
K 10000 virus /m? Assumed 





TABLE 4. Sensitivity Indices for Parameters of the model 2.1 





Parameter Index value Parameter Index value 





Ôp 


a d$ 


+0.00345 
+0.01530 
+0.33654 
+0.00811 
—0.30751 
—0.62107 
—0.58770 





v 
d 


Qe 
" 
Q 
Ho 


+0.38706 
+0.58770 
+0.25116 
—0.07953 
+0.00647 
—0.58770 
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The positive indices shows that the basic reproductive number increases as the value of the 
parameter increase. In Table 4, parameters a, b, o, Ac, Y, ^y, @ and d have positive indices. 
These parameters have a positive influence on the basic reproductive number, Ro. The 
contact rate between the susceptible population of village chicken and wild birds with con- 
taminated environment (d) is the most positive sensitive parameter of the model followed by 
the transmission rate between chronically infected village chicken and the susceptible chicken 
(v). The least positive sensitive index values is the transmission rate between chronically 
infected wild birds and the susceptible wild birds (9). The negative index values implies the 
inverse relationship between the basic reproduction number and the value of the parameter. 
Parameters ôe, ôb, K, u, and fl, have negative index values. According to these results, the 
Newcastle disease induced death rate (6,) in village chicken is the most negative sensitive 
parameter in the model. This tells us that more deaths of the village chicken during the 
outbreak of the disease reduce interactions between the infected population and the at risk 
population and hence lower the basic reproduction number. We observe also that the clear- 
ance rate (j,,) have the inverse relations with the Ro and thus increasing clearance of the 
NDV from the environments lowers the value of the basic reproduction number and the vise 
versa. The natural birth rate (u) is the least negative sensitive parameter in the model. 
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FiGURE 2. Shows (a) the effects of the clearance rate of NDV virus from 
the environment on Ro, (b) the effects of the shedding rate of NDV into the 
environment on Ro, (c) the effects of the contact rate of carrier wild birds and 
susceptible village chicken on Ro and (d) the effects of the contact rate between 
the susceptible population of village chicken and wild birds with contaminated 
environment on Ro. 


According to Figure 2 (a), increasing the clearance of the active Newcastle disease virus from 
the environment leads to the reduction of the value of Ro and hence reduces the rate of 
spread of the NDV in the avian population. Figure 2 (b) shows the direct relationship of the 
shedding rate of the NDV in the environment by the hosts on Ro. The more the shedding 
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the higher the degree of the NDV transmission in the population. Figure 2(c) shows the 
direct effects of the contact rate of the village chicken and wild birds with the contaminated 
environment with the NDV. More contacts during the epidemic leads to the increase in 
the rate of spread of the Newcastle disease. Figure 2(d) shows the direct relationship of 
the contact rate between the susceptible population of village chicken and wild birds with 
contaminated environment on Ro. Therefore in order to control the spread of Newcastle 
disease, the most positive sensitive parameters needs to be controlled. 


3. NUMERICAL SIMULATION OF THE BASIC MODEL 
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FIGURE 3. The dynamics of Newcastle disease in village chicken without any 
control. (a) shows susceptible population of the chicken (b) and (c) shows the 
exposed and infected population of the chicken and (d) shows the NDV in the 
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FIGURE 4. The dynamics of Newcastle disease in the wild birds. (a) shows 
susceptible population of the wild birds, (b) shows the exposed population of 
wild birds, (c) shows the infected population of the wild birds and (d) shows 
the carrier population of the wild birds 
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FIGURE 5. The transmission dynamics of Newcastle disease in village chick- 


en population, wild birds population and the NDV in the Environment dur- 
ing the outbreak in the absence of any control with initial state condi- 
tions Se(0) = 2000, Æ. (0) = 1200, 7e (0) = 1800, S, (0) = 3000, Æ (0) = 
1000, 7, (0) = 900, J, (0) = 700, H (0) = 5000 and the estimated parameter 
values % = 0.0083,a = 0.01,y = 0.067, u = 0.000548, 6. = 0.01984, ð, = 
0.025.d = 0.025, p = 0.9, = 0.001667, 155 = 0.00219, o = 0.02,b = 0,21 and 
Ro ~ 6.867. 


The Figure 3(a-d) and figure 4(a-d) show simulations of the dynamics of Newcastle disease 
in the village chicken and wild birds population in the endemic situation respectively. In 
figure 3(a), the susceptible village chicken population is exponentially decreasing due to the 
presence of ND in the population. Figure 3(b) represents the exposed population of the 
village chicken. The population decreases exponentially right from the first day after the 
infections. Normally infections takes place between day one to day six before the observable 
clinical signs of the disease on the individual chicken. Therefore as time goes on some chicken 
die naturally whilst others progress to the severe sickness class as shown in Figure 3(c). The 
population of the infected chicken increase for a short period of time due to new chicken 
which progress from the exposed class. The population of this class suddenly decreases 
exponentially due to the higher mortality rate of the ND. Figure 3(d) represents the NDV in 
the environment. We observe that during the first 100 days since infection the concentration 
of NDV in the environment increases to the carrying capacity due to higher shading rate of 
the NDV by the hosts. At its maximal point depending on the environmental conditions, 
the concentration of active virus decreases and thus lowers the transmissions of the disease. 
Figure 4(a) represents the population of susceptible wild birds in which we observe a sharp 
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decrease due to higher rate of infection of the NDV. Due to short window period of the 
ND, we observe a marginal increase of the population of exposed wild birds during the first 
few days of the disease outbreak which is followed by an exponential decline to the lowest 
level (see Figure 4(b)). The similar pattern is observed in Figure 4(c), as the population 
of the infected wild birds increases during the first few days since the disease outbreak and 
thereafter followed by an exponential decline to the lowest level due to the fatality of the 
ND. Figure 4(d) is a population of the survivors of the disease. This includes wild birds with 
immunity against ND in such a way that they remain carrier of the virus throughout the 
outbreak. Figure 5 depicts the dynamics of the three population (village chicken, wild birds 
and NDV in the environment) at the endemic situation when plotted in the same axis. 


4. DISCUSSION AND CONCLUSION 


In this paper, we developed and analyzed a mathematical model to elucidate the transmis- 
sion dynamics of Newcastle disease in village chicken which include vital components of the 
wild birds and the contaminated environment. The invariant region and positivity analysis 
of the model is presented which showed that the model is biologically meaningful and well 
posed in the feasible region. The basic reproduction number, Ro is computed and discussed. 
The sensitivity analysis results showed that, the shedding rate of NDV from the chronically 
infected village chicken, the infected wild birds and the carrier population of wild birds to 
the environment and the contact rate between susceptible chicken and wild birds with the 
environment (d), have direct impacts on the basic reproduction number, Ro. The increase in 
values for these parameters leads to the increase of the basic reproductive number, Ro whilst 
their decrease lowers the value of Rp. The clearance rate of the virus from the environment 
have also a great influence in the dynamics of the Newcastle disease among village chicken 
population. Increasing of the clearance rate, lowers the value of Ro. Hence, to control the 
rate of spread of Newcastle disease in the Village chicken it is of a paramount importance 
to limit the contact between wild birds and chicken and to increase the clearance rate of 
NDV in the environment through sanitation. However, more improvements to include the 
age structure and seasonal variations on the model are needed for further understanding of 
the dynamics of Newcastle disease in the village chicken population. Our model therefore, 
provides a baseline for the investigation of the control strategies of Newcastle disease among 
village chicken under a free range management system. 
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